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Abstract 

The fluctuations in the inflaton field at the end of inflation which seed the density perturbations 
are prepared in a pure quantum state. It is generally assumed that some physics causes this pure 
state to decohere so that it should be treated probabilistically. We show that the entanglement 
entropy between the universe inside our observable horizon and that outside our horizon is sufficient 
to do this. For the modes which are super-Hubble at the end of inflation, this entanglement entropy 
grows with volume inside the horizon, rather than with the horizon's area, and is proportional to 
the number of e-folds since Hubble crossing. 



1 



I. INTRODUCTION 



Inflation provides an elegant description for the large scale appearance of our universe 
(very large, approximately homogeneous and isothermic, and spatially flat to within mea- 
surement errors) [l|, 0, 0]. It also provides a mechanism to explain how the density perturba- 
tions which seeded structure formation (and therefore us) came to be (for a comprehensive 
review see IH). This occurs because, during inflation, the quantum state of (the approxi- 
mately harmonic oscillator /c-modes of) the inflaton field become highly squeezed 0. 

This squeezing, present on super-Hubble scales at the end of inflation, means that the 
inflaton is effectively taking on different values at widely separated points in space, which 
leads to an inhomogeneity in the temperature after reheating (or certain nonzero metric 
elements in a different gauge choice). 

One puzzle in this picture is that the state of the inflaton field is a quantum state. It 
is not really correct to say that it is "ahead" at some points and "behind" in others; its 
quantum state contains amplitudes to be ahead and behind at different points. Usually 
in treating the subsequent evolution of the system one replaces this quantum state with 
a classical probability distribution or ensemble with the same statistical properties. This 
artificially introduces a large entropy into the state, since a pure quantum state is being 
replaced with a statistical ensemble. It is an interesting question whether this procedure is 
really legitimate, and if so how it is justified. 

This issue has been mostly resolved by Polarski et al, who showed |6( that the differ- 
ence between performing calculations using the quantum squeezed states and the classical 
probabilistic distribution is exponentially suppressed in the number of e-folds of inflation. 
However it is still an interesting question whether the statistical treatment is in some way 
better justified than this. Are we actually compelled to consider the primordial fluctuations 
probabilistically? 

Burgess et al have argued that decoherence of the quantum state actually occurs due to 
interactions between the long-wavelength modes of the inflaton fields and shorter wavelength 
modes 0]. Assuming that these shorter wavelength modes are not observed, the description 
of the long wavelength modes should be made by integrating them out, resulting in an 
entanglement entropy. The calculation depends on the details of the interactions between 
long and short wavelength modes and on the non-observation of the short wavelength modes. 

We believe that it is unnecessary to consider this division of modes into long and short 
wavelength to show that the universe we observe will possess a large entanglement entropy 
and must be described statistically. This is because the post-inflationary universe has a 
horizon. We can only make observations of our universe within our horizon patch. But 
inflation sets up correlations between our causal patch and the universe "outside" that 
patch. Indeed, the basis which diagonalizes the wave functional at the end of inflation is 
the momentum basis, which is maximally nonlocal. Since it is impossible to observe the 
universe outside our horizon, we are obliged to integrate it out. The goal of this paper is to 
investigate the degree of entanglement with the universe outside our horizon. In particular 
we will compute the entropy of entanglement. 

Our method is based on that used by Srednicki [8|, who considered the entanglement 
entropy associated with cutting a spherical region out of space, as a toy model of a black 
hole's entropy. We extend his treatment to highly squeezed quantum states. Our main 
conclusion is that, for the highly squeezed case, the entanglement entropy grows with the 
volume of the visible universe, rather than with the surface of the boundary between the 



2 



observed and unobserved regions, as occurs for vacuum fluctuations. The entropy per mode 
scales with the log of the conformal time, which is the same as the number of e-foldings of 
inflation. This saturates a bound placed by Kiefer et al [9] on the entropy per mode which 
is allowed without destroying the phase information in the fluctuations, responsible for the 
acoustic peaks. Technically this conclusion applies only for the very longest wavelengths; 
if inflation took n e-folds to stretch our observed universe out of one Hubble patch, then 
the volume behavior we find is valid for fluctuations with wavelength more than 1/n of the 
current horizon length. (This volume dependence is consistent with the entropy conjecture 
of Bousso 13] because it applies only to the few infrared modes which were super-Hubble 
at the end of inflation and so became highly squeezed. The horizon area behavior argued 
by Bousso is dominated by the much more numerous shorter wavelength modes.) 

The entropy per mode which we compute is an inevitable result of the horizon structure 
of the post- inflationary universe and does not rely on any assumptions about the interactions 
of the inflaton field with other (short- wavelength) degrees of freedom. It is therefore robust, 
and sufficient to explain the statistical nature of the universe we observe. 



II. SQUEEZED STATES 

A squeezed harmonic oscillator state is generated from the vacuum state by a squeezing 
operator, which is a function of two variables, the amplitude of squeezing r and the phase 
$. The squeezing operator for a momentum-space mode of wave vector k of a quantum field 
is generically by Eq. below. 

During inflation, the state of a noninteracting, light scalar field can be described by the 
quantum state of each momentum-space state. The cosmological expansion causes the mo- 
mentum space states with wavelengths longer than the Hubble scale become highly squeezed, 
in a manner which has been worked out by Albrecht et al |5]. They find that the evolution 
equations for r k and in terms of proper time 77 are 
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Here U is the time evolution operator of the inflaton and R(Q k ) is an overall phase. Note 
that these solutions are for the exponentially expanding deSitter stage, in the Bunch Davies 



vacuum. 
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Evaluating the state of this fc-mode in the field basis gives [11] 

= (<f>(k,v)A(-k,r])\S\0) 

= iV fe exp(-|<XM)| 2 ^) (7) 
1 +isin(2$ife)sinh(2r fc ) 
cosh(2rfc) — cos(2$fc) sinh(2rfc) 

Here (j)(k,rj) is the scalar field, called y(k,rj) in is a normalizing constant; we will 

not keep track of such normalizations, ensuring the correct normalization at the end by 
scaling the density matrix to have trace 1. The fact that the wavefunction can be put in 
this simple, Gaussian form means that we can adapt the approach Ref. (s| used to calculate 
the entropy of a sphere of space to calculating the entropy of the squeezed states. 



III. INTEGRATING OUT POSITION-SPACE MODES 

As discussed, we want to compute the entropy of entanglement of the visible part of 
the universe (that which lies within our horizon) when the invisible part (space outside our 
horizon) is integrated out. We are only interested in the behavior and entropy of relatively 
long wavelength modes, those responsible for the density perturbations on scales which are 
linear today or can be probed at a time scale when they were linear. At the end of inflation 
and in comoving coordinates, these modes had wavelengths much larger than the Hubble 
scale and so were all highly squeezed. To describe only these long wavelength modes, we 
should implement some infrared regulator, which for convenience we take to be a lattice 
discretization of space. (We do not think that the choice of regulator is important to 
our results.) We will also make our system finite by imposing boundaries far outside the 
observable horizon. We will seek the limit where the boundaries are moved to infinity and 
will study the dependence on the infrared regulation. As a further simplification we consider 
a 3-dimensional system where the visible universe is a cube rather than a sphere; again the 
qualitative behavior should be the same as for the spherical case. As a test case we also 
consider a 1-dimensional system. 

This effectively reduces the problem to a quantum mechanical problem of the product 
of a large number of harmonic oscillators. We will find the wave function and from it the 
density matrix, which is most easily done in the momentum basis. Then we will transform 
to the position basis and integrate over the coordinates outside the horizon. The procedure 
is similar to that of Srednicki [8| , except that we consider highly squeezed SHO states, rather 
than vacuum states. 

In one dimension, the discretized box Hamiltonian for a free scalar field is 



H 



N 

E 



2<j)( Xj ) + 4>{Xj-l)\ + 



(9) 



where I is a lattice spacing, Nl is the extent of the universe, and 4>(xj 
N generalized coordinates and conjugate momenta. Here (<f>j+i — 2<f)j + 4>j-i) is a discrete 



and TTj are the 
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approximation of / 2 V 2 0, and we switch to matrix notation in the last expression. We take 
0(0) = = 0(iV+l), which enforces Dirichlet boundary conditions. This choice of boundary 
conditions should be immaterial if we take an appropriate large volume limit. We have tested 
this by using mixed Dirichlet-Neuman and antiperiodic boundary conditions. 

The matrix M is diagonalized by a sine transform. Explicitly, defining the matrix U as 



iV + 1 \N+l 
the matrix M becomes M = Wu 2 U, with 



u = Diag[Wj] , ^^(^L^j. (11) 

The wave function for the scalar field is easiest to express in this momentum basis. 
Defining <pj{k) = U^cp^x), then up to a normalization constant, the squeezed state in the k 
space field basis is 

N 

(<p\s) = n(<f> k \s k ) (12) 



fc=i 

N 



J^[exp(-fi 



k=l 



= exp(-0(fc) t fi0(A;)), 

where <f>(k) is a vector of N k modes and Q is the diagonal matrix with elements flk given 
by Eq. (jSJ) with u>k replacing k. 

We can express this in the position basis using the unitary transformation matrix U found 
earlier. Defining M = WQU, the state in position basis is 

((j)\S) = exp -<j)(x)W(f)(x) , (13) 

giving the density matrix 

p(0(x),0(x)') = (0|5)(5|0') 

= exp - (<j)(x)W<f)(x) + (P(x)'W^(x)') . (14) 



With the density matrix in this form it is now possible to calculate the partial trace of the 
physical modes outside the visible universe using an approach similar to that of However 
before continuing it is important to note that the method used there must be altered slightly 
because M is not Hermitian, that is, Q is not real. 

Now we subdivide the iV lattice points into m points making up the visible universe and 
n = N — m modes lying outside our horizon, which we therefore want to integrate out. In 
the 1-D case the m modes are taken to be bounded on each side by | modes, and in the 
3-D case the visible universe is taken to be the central modes of an iV x iV x iV cube. 

We should therefore write the matrices U and Q in block form as follows: 



u, u 2 



L rn 
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From here it is straightforward to see that 



[y ] % ] ] 



A R + iA c C R + iC c 
C R + iCt B R + iB c 





y 




X 



(17) 



Where 



A R 


= U{ Re [ttnpx + U\ Re [n m ]U 4 


(18) 


A c 


= U\ Im [n n ]U x + U\ Im [Q m ]U 4 


(19) 


C R 


= U\ Re [Q n ]U 4 + U\ Re [Q m }U 2 


(20) 


c c 


= U\ Im [Q n ]U 4 + U\ Im [n m )u 2 


(21) 


Br 


= UlRe[n n )U 3 + U^Re[n m )U 2 


(22) 


B c 


= U\ Im [n n )u 3 + U\ Im [n m )u 2 


(23) 


B 


= B R + iB c . 


(24) 



Note that A RtC , B Rc are all Hermitian matrices and we have split the N modes into a y 
vector of n modes and an x vector of m modes. We will be tracing over the y modes. The 
partial trace of a density matrix can be found by integrating over the appropriate modes; 
that is, 



Pproj {x, X ) 



dy p(x,y;x',y). 



(25) 



Replacing y' with y, the density matrix to be integrated is 



p = exp 



2y ] A R y + y ] [C R {x+x') + iC c {x-x')\ + \{x+x') ] C R + i{x-x') ] C ] c ]y 
+x^Bx + x'^B^x' 



(26) 



The Hermitian matrix A R can be written as a unitary transformation of a diagonal matrix, 
A R = V^XV. After completing the square for y and carrying out the Gaussian y integration, 
the density matrix becomes 

p proj = exp - (x\K - J - i(G + G f ) + B)x + x'\K - J + i(G + G ] ) + B ] )x ! 
-x\K + J + i{G- Gfi))x' - x'\K +J-i(G- G f ))x' 



= exp — ^Tx + x ^T^x' — x^5\x' 



x ^5 2 x 



where 



A" 



G = ClV ] VC 

2 

K = [VC C + (VC t 

J ee C R V^VC J 

T = K-J-i{G + G^) 
5 1 = K + J + i(G-G^) 
8 2 = K + J-i{G- G ] ) 

In the case where Q is purely real, this reduces to the result of |8| 



R 
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(28) 
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(31) 
(32) 
(33) 
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IV. ENTROPY, EIGENVALUES, AND EIGENVECTORS 



A. Diagonalization of p into normal modes 

We now have the density matrix for the visible part of the universe. The question is, can 
this be represented by a stochastic distribution, or is its quantum nature still observable? To 
determine this we calculate the amount of decoherence of the density matrix by computing 
its entropy. To do that we need to find the eigenvalues and if possible the eigenvectors of 
the density matrix. 

Since 5\ and 5 2 are both Hermitian, we can further simplify the mixed term in the density 
matrix: 

x^Six' + x' ] 5 2 x = z^Xsiz' + z' ] Q5 2 Q ] z (34) 
= z^Xsxz + z^Q6 2 Q^z = x'\5\ + 5 2 )x 

= 2x'\K + J)x, 

where Q^5\Q = A51 =diagonal, and x = Qz. 

Now we split r into its Hermitian and ant i- Hermitian parts: 

T = A + iD, (35) 
A = K-J + B R , (36) 
D = -G-G^ + B c . (37) 

This gives us a density matrix of the form: 

Pproj = exp — (x^Ax + X ,] Ax') + i(x ] Dx - x' ] Dx') - 2x\A - (B r - 2J))x' . (38) 

We know that A is positive definite since the density matrix is normalizable. Diagonalizing 
A = WAW^ and dividing by a/A gives 



Pproj — exp 



(v)u + u'V) + i(v)D'u - u ] D'u') - 2« f (l - A~^W\B r - 2J)W Ar^)u' 
where x = k l l 2 Wu. Diagonalizing A~^W^(B R — 2J)Wh~^ gives 

Pproj = ( fl i/^exp(-[(«J + vf) - 2^(1 - e>;])] eM-^D"v - v'^D"v')) 



(39) 



\3=1 



Ylpjivj- v'j) exp(-i(vW"v - v'^D"v')) , (40) 

where v is in the new diagonal basis, and €j are the eigenvalues of h~^W\B R — 2J)H^A~^. 
The second equation defines pj(vj] v'j) = \ ^ exp(— [(i>| + v'f) — 2vj(l — ej)v'j]). 



B. Proof of unimportance of complex phase 

Besides the D" term, this density matrix in Eq. fT40|) is in the same form as found in 
We now show that the D" term does not change the eigenvalues and only modifies the 
eigenvectors in a trivial way. 
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Write the eigenfunctions of p{vi,v'-)i as f n (i);i, satisfying 

Pi(Vi, V' i )f n (i);i(v' i )dv' i = Pn(i);ifn(t);i(Vi) , (41) 



with p n (i);i the associated eigenvalue. Then the eigenvectors and eigenvalues of our density 
matrix are a D dependent phase factor times products of these, namely 

m 

9n { i )M 2),.Mm)(v) = e { ~ ivtD " v) J] Li J, AO) , (42) 

3=1 

m 

Pn(l),n(2),..,n(m) = IlPn(j)a ' ( 43 ) 
3=1 

This is easy to check: 



/oo poo I '"- \ "> 

dv> Ppmi (v,v>)g(v>) = / dv> mpitort) e H(^''-^V)) e (- i „^v)Jj /) 

/oo / m 

OO V „-_1 



3=1 



= exp(-iv^D"v) lYlPn(j)fn(j)A 

\j=l 

= Pn(l),n(2),..,n(m)9{v) ■ (44) 

Therefore the complex phase does not actually affect the eigenvalues of the density matrix, 
and only adds an overall phase to the eigenvectors. However this complex phase will be 
important when taking further partial traces. 



C. Entropy per mode 

It was shown in [8J] that the density matrix pj(vj;v'j) found in Eq. (f40|) has eigenvectors 
and eigenvalues of 

fn(vj) = = H n (a-*Vj) exp (-^vf) , (45) 

a/ an\2 n y/7T v ' \ l J 

Pn = ( 1-Or , (46) 

a = V2e - e 2 , (47) 

< - rrf ■ (48 » 

In other words, up to the phase factor which we just showed is irrelevant, there exists a 
basis in which the density matrix behaves like the product of density matrices for individual 
harmonic oscillators, with Poisson distributed occupancies of states which are the same as 
the conventional harmonic oscillator number states. 
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The entropy is the sum of the entropies for the normal mode density matrices. The 
entropy of the Poisson distribution described by p n = (1 — is 

^ = _l 0g (l_£)_ J_i og (£). (49) 

As we will see, we will be interested in the limit where €j <C 1, in which case £ ~ 1 — v2e 
and the entropy is 

s = E%^E( 1 -^ lo g( 2 ^)) • ( 5 °) 



D. Large Squeezing Limit 



For large squeezing, ~ k{k 2 rf + ik\i]\), which is small and dominated by the imaginary 
part. The overall scale of fl k is immaterial; multiplying all Q k by a common factor scales 
out in determining the 6j and therefore the entropy. However, the k dependence and the 
relative size of the real and imaginary parts is relevant. Since krj is exponentially small in 
the number of e-foldings, we should study the scaling of the eigenvalues ej in the small k\r)\ 
limit. 

Ultimately the entropy depends on the sizes of the eigenvalues 6j of the matrix 
A~5W'(Br — 2J)WA~2 in Eq. ()40l) . Therefore we should study the scaling behavior in 
rj of the pieces of this matrix. It is convenient to multiply Qk by l/rj, with I the lattice 
spacing (shortest wavelength under consideration), so 

tt k ( V )^k(^ik) 2 (^+i(ik)y (51) 

In the small ^ liinit , provided that the eigenvalues of K are all nonzero, one may neglect 
Br — J in Eq. (1361) . In this case the rj dependence scales simply; the matrix A~zW\Br — 

i i 1 2 

2J)WA~2, and therefore the £j, scale as ? . Defining e(i)j as the set of eigenvalues arising 
by setting ^ = 1 and setting A = K in Eq. (|36|) . the approximate values of the ej are 



T] 2 



(1),. (52) 



Therefore the entropy is approximately 

5^^1n^- + l-iln2-iln6(l),. (53) 



\V\ 



This demonstrates that, for rj/l parametrically small (modes which were well outside the 
horizon at the end of inflation), the entropy is extensive in the number of modes and scales 
as the number of e-folds by which the mode was outside the horizon (the log of the level of 
squeezing) . 

Note however that this conclusion relies on the e(l)j not being exponentially large. Some 
of these eigenvalues could potentially be exponentially large if there is another parametrically 
large number in the game, such as the number m of points inside the horizon. Therefore 
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our conclusion that the entropy scales as the volume and as — In | krj | is only robust if we 
consider fluctuations on very long scales, m < — In \rj/l\. 

We have been unable to find a closed-form solution for the eigenvalues e(l)j. Instead we 
have performed numerical experiments to check how the residual entropy in Eq. (1531) scales 
with volume. Recall that the choice of the number of points inside the visible universe 
corresponds to a selection of what scales of perturbations to consider, while we should seek 
the limit that the volume outside the visible universe becomes large. 

We studied the entropy in both the one dimensional and three dimensional cases. The 
matrix manipulations rapidly become numerically unstable and were only able to achieve 
a large value of m for the 1-dimensional case. For small m (an analysis only considering 
modes of near-horizon wavelength), we find as expected that all eigenvalues e(l)j are 0(1) 
and therefore the entropy grows with the number of modes and scales with the number 
of e-folds. When we increase m, we observe that the mean value ^ Sj=i m me (l).j grows 
roughly linearly with m; the largest eigenvector Vj is one which avoids the boundary. This 
suggests (though we have not been able to show it analytically) that the volume behavior 
we obtain will break down when one considers wavelengths shorter than the current horizon 
by more than the number of e-foldings of inflation. 

Return to the longest wavelength modes, where we have shown that the entropy per 
mode is ~ Refolds- This is exactly what Kiefer et al [qJ predicted to be the upper bound 
on the entropy of the modes. They argue that if the entropy were any higher, then the 
experimentally observed acoustic oscillations of the CMB fluctuations would not be seen. 
Therefore, so far as the modes which can actually be observed are concerned, the decoherence 
arising from entanglement with the (unobservable) region outside our horizon is sufficient 
to force a probabilistic description of the density fluctuations, regardless of whether or not 
there are additional mechanisms for decoherence. 

V. CONCLUSION 

The post-inflationary causal structure of the Universe has a horizon, meaning that we 
can only observe physics in some finite domain of the universe. However, inflation induces a 
large entanglement between the IR field modes in this region and those outside our horizon. 
We have calculated the entanglement entropy for the IR modes responsible for density 
perturbation seeds. The entropy scales with the number of modes in the visible universe, 
that is, with the volume of the universe times the cube of the fc-vector at which density 
perturbations are measured. The entropy per mode obeys the simple rule s ~ r, where r is 
the squeezing factor (the number of e-foldings of inflation during which the relevant mode 
was outside the Hubble radius, of order 60). This far exceeds the cosmic variance limit 
on the amount of information available in the density perturbations. In fact, it saturates 
the upper bound predicted by Kiefer et al jsj. This result holds for the longest wavelength 
modes, those with HqX > 1/r, with Hq the current Hubble constant (roughly the inverse of 
the horizon length). 

This entanglement entropy is large enough to render unmeasurable any quantum coher- 
ence in the initial squeezed state set up by inflation and to justify a probabilistic description 
of the initial seeds of the density fluctuations. 
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